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ABSTRACT 


A COMPARISON OP SIMULATED AND EXPERIMENTAL WAVE SPECTRA 
IN THE NEARSHORE REGION 

W. Douglas Morris 
Old Dominion University, 1979 
Director: Dr. Chester E. Grosch 

The increasing use of the world T s continental shelf and shoreline 
for resource > recreation, and commercial activity will require im- 
proved abilities to monitor and predict the wave climate of these areas 
for safe and efficient mapagemqnt of the nearshpre region. For predic- 
tive models to be used in this process,, they must first be verified by 
comparison with in situ wave measurements:. This paper compares £116 wave 
energy spectra from one suoh model, with spectra from continuous wave 
height measurements of a wavefield mqving from deepwater to the shore- 
line after passage of tropical storm Amy on July 2, 1975* Discussion 
of the theoretical model, data gathering techniques, historical develop- 
ment and the analysis procedures are also included. 

The results of comparing the experimental spectra measured over 
the shelf with spectra measured in deepwater conditions indicated that 
statistically significant changes occurred over the shelf region that 
could be attributed to shoaling and refraction effects. The comparison 
of the simulated spectra with the experimental agreed well at the 80 % 
confidence level for all but the most inshore data segment. The total 
energy also showed similar agreement*, however, the theoretical frequency 

/ 

ORIGINAL IS 

OF POOR QUALITY 



at which peak spectral density occurred did not follow the decreasing 
trend seen in the experimental data, pui delines are proposed for 
determining the type wavefield to which the modeling technique can 
test he applied, and caution is given pn the selection of the digiti- 
zation rate necessary to ensure sufficient detail for analysis of 
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INTRODUCTION 


The increasing use of satellites for synoptic monitoring of the 
ocean surface conditions will furnish information that can be used to 
better plan and utilize the nearshore regions. A combination of satel- 
lite sensors will measure wave height , direction , wind speed and provide 
data that can be used for the computation of wave spectra (McCandless, 
197^)* However, these measurements will represent a broad-scale look 
at the ocean surface in offshore areas where the surface conditions 
are relatively homogeneous over the swath coverage of the satellite. 

i 

By using this data obtained from remote sensing in conjunction with 
theoretical computer modeling techniques, it should be possible to ex- 
tend the useful range of the satellite data further into the coastal 
zone and to provide more detailed short-term predictions of the sea 
conditions in the nearshore region. One such technique requires only 
the initial deep-water wave height, direction and wave period for use 
in a linear wave refraction model to determine changes in wave charac- 
teristics as they propogate across the shelf. Such capabilities, when 
combined with satellite data would make possible an extended use of the 
satellite data in coastal zone planning and management , for such as 
coastal disaster warning, coastal protection and development, and off- 
shore siting for facilities. Before such a modeling technique could be 
used as an operational management tool, the results of the model would 
have to be verified by comparison with measured ocean wave data. 

The purpose of this research was to compare simulated wave energy 
spectra, with wave energy spectra from experimental data. The experimental 
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data were recorded off the coast of North Carolina on July 2, 1975, 
after passage of tropical storm Amy. Wave energy spectra were computed 
from the experimental data using the .method of Blackman and Tukey (1958). 
The simulated spectra were estimated using a theoretical linear refrac- 
tion computer model and a technique based on the method of Pierson, 

•4 

A 

Neumann, and James (1955) for constructing wave spectra. 
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HISTORY MD RELATED WORK 

As waves cross from deep water into shallow water their speed and 
direction change as an interaction occurs between the wave form and the 
sea floor. That portion of the wave which first encounters shallow 
water is slowed , while that portion that remains in deep water continues 
at a speed greater than the shallow water wave speed. This causes a 
bending and reshaping of the wave crest known as refraction. These 
changing wave patterns affect the distribution of wave energy over the 
shelf and along the shoreline. 

Theoretical Models 

Prediction models have been developed which, given deep water 
input parameters, can be used to determine these changes in wave energy. 
Early models such as that developed by Johnson and O'Brien (19^6), used 
graphical construction of each wave crest as it advanced inshore. Later, 
this model was expanded to construct wave ray diagrams without drawing 
the wave crest (Johnson, O'Brien, and Issacs, 1948). This was still a 
graphical technique with the wave rays perpendicular to the wave crest. 
Graphical processes like these were later incorporated, into computer 
programs by Griswold (1963) S Mehr (1962)-, and Harrison and Wilson ( 196 ^), 
but these programs used a spatial grid of wave speeds to calculate the 
path of the wave ray and thus limited its usefulness by requiring a 
new grid for each new wave' period. A newer model developed by Wilson 
( 1966 ) , eliminated many of the problems associated with the wave speed 
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grid "by using a grid of water depths in computing wave ray position. 
Dobson ( 1967 ) > combined wave ray construction using a grid of water 
depths with a finite-difference solution to the equations of wave in- 
tensity as developed by Munk and Arthur (1952), to compute wave height. 

With the advent of the computerized refraction model, it became 
practical to compute the effects of refraction on ocean waves in the 
manner first proposed by Pierson, Newman n, and James (1955)- Their 
theory called for the refraction of a large number of simple sine waves 
from various directions and with various periods about some point of 
interest in shallow water. From these results, plots could be construc- 
ted of two functions which show the variations in refraction and direc- 
tion of waves about this location. The spectrum of the deepwater waves 
that would be refracted would then be divided into an approximating sum 
of simple sine waves, with a frequency, direction and a fraction of the 
total energy of— the spectrum characterizing each sine wave. These values 
would then be used with the previously constructed refraction and direc- 
tion charts to compute the incremental energy and direction at the point 
of interest for each deepwater sine wave. By then summing the results 
of each sine wave from the deepwater spectrum both the wave height and 
direction at the shallow water location could be determined. This is 
possible because theoretical considerations indicate that, even for 
spectral frequencies separated by only a small amount, the individual 
sine waves behave independently of each other as far as refraction is 
concerned. With modifications, this principle has been applied to 
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several recent papers. Chao (197*0 used this technique to predict 
the refracted wave spectrum near the Chesapeake Bay entrance, and 
Poole (1976b) used this technique to compare remotely sensed wave spec- 
tra, with that computed using a wave refraction model similar to 
Dobson 1 s . 


Instrumentation 

For a number of years oceanographers have been searching for ways 
to improve recordings of wave height distributions and a number of in 
situ instruments have been developed, but it was not until the early 
fifties that the first wave measurements were made from an aircraft 
using a wave profiling radar (Ross, Peloquin, and Shell, 1968). Follow- 
ing the advent of lasers, the laser wave follower profilometer was 
developed in the late sixties (Ross, Peloquin and Shiel, 1968). The 
profilometer has since been used in other studies by Schule, Simpson 
and DeLeonibus ( 1971 ) and Poole, ( 1976 b). 

'Experimental Analysis 

Because- of the random nature of ocean waves , a statistical approach 
is desirable in defining the properties of the sea surface! It was in 
1952 that Pierson (See Kinsman, 1965), building on the theories of Tukey 
and Hamming (See Kinsman, 1965 ) on power spectra in the field of elec- 
tronics, applied these theories to the study of ocean waves. Newmann 
(See Kinsman, 1965) applied the power spectrum theories to ocean wave 
time histories to obtain the wave energy spectrum as a function of mean 
wind. St. Denis and Pierson ( 1953 ) first addressed the problem of using 
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the spectrum of a wave record made from a moving platform and mapping 
this into its counterpart in a stationary* reference frame. This 
spectrum is valid only in the moving reference frame, and thus gives a 
distorted estimate of what would "be the true frequency estimate in the 
stationary reference frame. Cartwright (1963) used the theories dis- 
cussed in St. Denis and Pierson to map the spectrum from the moving 
reference frame into its counterpart in a stationary reference frame. 

These transformation relationships were used in a study on wave genera- 
tion by Barnett and Wilkenson (1967), Ross, Peloquen, and Shiel (1968) 
and Schule, Simpson, and De Leonibus ( 1971 ) • These latter two reports 
used a laser profilometer to measure wave height histories where the 
first report used radar. 

Poole (1976a) showed that the deep water approximation to the 
linear wave dispersion relationship, as had been previously used, intro- 
duced errors when applied in areas of shallow or intermediate water depths 
where refraction takes place. Poole developed transformation equations 
valid for all depth zones and used these in the study of a wave field 
(Poole, 1976 h) where a comparison was made between remotely sensed wave 
spectra and spectra computed by using a wave refraction computer model. 
However, most of the energy of the wave field was at frequencies higher 
than the frequencies for which refraction was predicted. In addition, 
the variations among the remotely sensed spectra were statistically in- 
significant and thus the variation could not be attributed to any physical 


causes . 
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The present study was of a higher energy wave field associated with 
the passage of a tropical storm. The frequencies at which the energy 
would he concentrated could he expected to he those for which refraction 
effects would he predicted. In addition, longer wave lengths and periods 
associated with storm conditions would produce stronger shoaling and 
refraction effects which increase the likelihood of significant changes 
in the levels of the experimental spectra. For these reasons, the 
methods employed in Poole (I976h) were applied to the wave data from 
tropical storm Amy. 



8 


GOVERNING EQUATIONS 

Wave energy spectra is a measure of the energy distribution with 
frequency. Changes in the spectral shape are indicative of a shift- 
ing of the wave energies to different frequency bands . As waves cross 
the shelf region the change in spectral shape can be used to compare the 
redistribution of the wave energy. The spectra can be computed either 
from measured wave heights or from the results of a wave refraction 
model. The equations governing the computation of energy spectra are 
presented as follows: 

Spectra from Measured Data 

The energy in the ocean waves is a function of the fluctuations in 
water height. This fluctuation is a random process and as such is 
describable by the laws of probability. In order to generate a wave 
energy spectrum from records of finite length, it is necessary to assume 
that the process is both stationary and ergodic (Kinsman, 1965)* Using 
a time history of the surface elevation variation about the mean, a wave 
spectrum can be computed by first taking the autocorrelation of the data 
record. The autocorrelation function is a measure of the degree of 
parallelism between a given time series n { t.) and the same series at a 
displacement x. From Bath ( 197*0 , the autocorrelation is computed by 

N-x 

C(t) = ~ 2 U(t)r](t+T) (1) 

t=l 


T = 0,1,2,. . .M 
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where M is the maximum number of lags- 

The raw spectral density estimates are then computed by taking the 
Fourier transform of the autocorrelation function using the following 
equation : 


C(t) COS -jjjr- + C(M) cos ilir 


which divides the spectra into M frequency bins of equal width (1/2 At). 
In transforming the autocorrelation function, these data have in effect 
been passed through a lag window (Bath* 197*0 and the undesirable spectral 
effects which are thus introduced are then compensated for by smoothing 
in the frequency domain. The smoothed spectral density estimates are 
obtained by applying the Hamming-Tukey smoothing equation to equation (2). 

= -23S(£ - 1) + •54S(A) + -23SU + 1) (3) 

These spectral density estimates represent the one-side cosine trans- 
form values. 

* In the case where the data are recorded in a moving reference 
frame, the frequency is shifted from that which would be seen in 
a fixed reference frame. In this case, the frequency associated with 
each spectral bin represents an apparent frequency as seen from the 
moving platform. In order to map the apparent spectrum into a stationary 
reference frame, the unidirectional transformation technique developed 
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by Poole (1976a) is used. This technique utilizes the fact that an 
equivalent amount of energy density must reside in both the apparent 
frequency bin and the transformed frequency bin. In the stationary 
reference frame, the width of the frequency bins vary, and the corres- 
ponding spectral densities must also change to maintain a constant 
amount of energy density. Thus the transformation is made in pairs, 
such that for each apparent frequency and spectral density pair in the 
moving reference frame, there will be a corresponding frequency and 
spectral density pair in the stationary reference frame. The transforma- 
tion uses specific values of aircraft velocity, apparent frequency, and 
averaged water depth for the region of study, to find a unique iterative 
solution for the wave number k (equation 4) under the assumption that 
the speed of the platform from which the measurements were made was 

kV - (gk tarih kd) 1 ^ 2 - w = 0 (k) 

a 

greater than the speed of the fastest wave. The corresponding frequency 
and spectral density can then be confuted in the stationary reference 
frame by using equations (5) and (6) 

(i) 2 = gk tanh kd (5) 


S(w) = S (a) ) 
a a 


2d) Y/g 

tanh. kd + kd seeh 2 kd 


( 6 ) 
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Simulated Spectra 

Wave energy spectra- can be constructed using wave refraction 
models as discussed in Pierson, Neumann, and James (1955 ) 5 Chao (197*0, 
and Poole (1976b). These refraction models predict the changing wave 
characteristics and crest patterns produced by waves as they interact 
with the bottom topography (Goldsmith, et al. 197*0* As the waves move 
from deep water into shallower water, the wave speed is slowed and the 
wave shape becomes shorter and steeper. Since the process Is not uni- 
form, but varies with the changing depths, the varying wave speed along 
the crest causes the bending of the wave front* This process is analogous 
to the refraction of light rays in optics and is governed by Snell ! s law: 


sin a 2 c 2 

sin a ± c^ 

where and are the angles 

between successive wave fronts and 
the corresponding bottom contours as 
shown in the accompanying sketch, 
and c^ and c^ are the respec- 
tive wave celerities associated with 
these wave fronts. Assuming steady 
state conditions and a progressive 
sinusoidal wave with small steepness 
and constant frequency, linear wave 


(7) 
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theory can he used to compute the local wave celerity from the dispersion 
relationship 


c 2 = tanb 
2ir 


2ird 

L 


(8) 


Using a deepwater direction and period, the local phase speed and 
wave curvature can he computed. These values are then used to compute 
a new phase speed and curvature for the next point at some fixed time 
step and hy this means the wave pattern is advanced shoreward. For a 
more detailed discussion of this technique see Dohson (1967), Goldsmith, 
et al. (197^) » and Poole, et al. (1977)* 

The result of refraction is to change wave height as well as 
wave direction. By assuming that no energy flows laterally along the 
wave crest , and that no energy is gained or lost due to reflection , 
percolation, or bottom friction, then the wave height can he determined 
from , 


H 

H 


K K 


r s 


where K is the refraction coefficient., 
r 



and K g is the shoaling coefficient 


(9) 


( 10 ) 
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K 


i 


2 cosh, kd 
2kd + sinh kd 


( 11 ) 


Since the wave energy is proportional to the square of the wave 
height, equation (9) can he used to compute a normalized energy ampli- 
fication function as follows: 


F(f) = l^y = K^(f ) K 2 s (f) . (12) 

o 

The amplification function F(f), is computed for the full frequency 
range corresponding to each spectral region using the refraction model 
to compute the refraction .and shoaling coefficients for the monochromatic 
wave frequencies of the reference spectrum. Combining these values with 
the spectral densities of the reference spectrum S Q (f), the inshore 
simulated spectra can then be determined from equation (13) as a 
function of wave frequency (Poole, 1976b). 


S(f) = S Q (f) F(f) 


(13) 



PROCEDURES 


Data Collection 

? 

The wave data used in this study are from wave measurements taken 

t 

on Jul^f 2, 19T5 j after .passage of tr'opical storm Any off the coast of 

i * * 

* P 

North Carolina. At, the time of the experiment (1600-1630 d. s.t.) 

V 

Amy’s storm center was approximately 512 n.mi. offshore and had maximum 
sustained winds of 4 60 knots. The ground track of the storm is presented 
in figure 1, and position data in table 1. During the U8 hour period 
prior to the data gathering flight* the storm moved further out to 
sea in a relatively constant direction from the test site. Figure 2 
is a photograph from the Synchronous Meteorological Satellite-1 ( SMS-1 ) 
showing the position and size of the storm near flight time. 

Data were taken as a part of a 
flight experiment using a NASA C-5^ 
aircraft stationed at Wallops Island* 

Virginia. The data leg (see sketch) 
began approximately 8l n.mi. off- 
shore at an altitude of 152 meters 
and flown in the onshore direction 
along the path of swell propaga- 
tion* as determined from pre flight 
wave reports and visual observation 
from the aircraft. A continuous 
surface profile was measured along 
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this line of flight using a laser wave profilometer and recorded on 
analog tape. Bata gathering was terminated at the shoreline over 
Jennette's Pier, Nags Head, North Carolina. An inertial navigation 
system was used to periodically record the flight position as presented 
in figure 3. Visibility was 7 miles or better in the experimental area 
with little cloudiness. Surface winds were from the Northwest at 5 to 
10 knots. 

The laser profilometer used to record the wave height time history 
was furnished and operated by the Naval Oceanographic Office (N00). The 
instrument was a Spectra Physics GeocLolite-3 ranging system utilizing 
a constant wavelength helium-neon laser with a wavelength of 632.8 nm, 
and had an output of 25 millivolts (Ross, Peloquin, and Shiel, 1968). 

Experimental Analysis 

The surface elevation time history was recorded and preprocessed 

k * 1 

r 1 * i * * 

by the WOO. Raw data were • measured by a laser profilometer and recorded 
as an analog signal. This record was then digitized at a 0.02 second 
interval and divided into 10 consecutive data segments. From the off- 
shore position* the first nine segments were each 3 minutes (9000 samples) 
long and the tenth was 1-727 minutes (5181 samples) long. The 10th 
segment was shorter as the data record was terminated at the beach. 

Thus 10 quasi -homogeneous data sets were available for analysis. Data 
were recorded in terms of voltage and converted to altitude measurements 
by using an WOO supplied data factor. A Dinette high pass filter 
(Linette, 1961) was used to remove low frequency (less than 0.15 hertz) 
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aircraft vertical motions from the data* Roll and pitch motions ve re 
assumed 1 to have a negligible effect on the data record when an 
aircraft of this size is used (Barnett and Wilkerson, 1967). 

For each data segment, the mean surface elevation was computed and 
subtracted from the data record. Subsequent computations were based on 
wave fluctions about the mean as it is the variability of the fluctua- 
tions which is the statistically stable quantity to be examined (iCinsman, 
1965)* The data record then represented surface elevation variations 
with time in the aircraft reference frame. This moving reference frame 
provided a compressed data record from that which would be observed 
in a stationary reference frame. This apparent wave spectrum was first 
computed and then transformed into the stationary reference frame to 
obtain the experimental spectrum. 

The apparent spectrum was computed using the autocorrelation 
function in equation (l) for each data segment. The number of lags M 
for this equation must be chosen as a compromise between the desired 
degree of correlation stability and spectral resolution. A large value 
of M increases the resolution but decreases the number of terms used 
to compute the correlation which might lead to fluctuations in the 
values of the autocorrelation functions. Blackman and Tukey (1958) 
recommend that M be no more than 5 to 10 percent of N for power 
spectrum computations, but aside from this restriction the value used 
becomes a matter of judgment. Figure b shows the autocorrelation as a 
function of T for 150 lags. For most spectra this represents ll 4 


degrees of freedom. No further significant changes in the autocorrelation 
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function was seen when larger values of M were tried. The Fourier 
cosine transform of this autocorrelation function was then used to 
compute the raw power spectrum according to equation ( 2 ) , and the 
values then smoothed by the Hamming-Tukey equation ( 3 ). Then these 
smoothed apparent spectral values were transformed to the stationary 
reference frame by using the unidirectional transformation technique in 
equations ( 4 ) , ( 5 ) and ( 6 ,). The’ apparent frequency values were 
calculated from the raw data, the aircraft velocity was 83.8 m/s and 
the average depth was determined by averaging the depth along the 
flight line for each segment. The resulting experimental wave spectra 
in the stationary reference frame are presented in figure 5 for the 
first 4 data segments. Since these first 4 spectra are each from 
measurements made under similar deepwater conditions, they may be 
regarded as statistically independent random variables and can be added 
together to improve the accuracy and statistical stability of the 
resultant spectrum (Lumley and Panofsky, 1964). This is then known 
as the deepwater reference spectrum and is presented in figure 6 . The 
remaining 6 inshore experimental spectra are presented in figure J. 

Theoretical Analysis 

A modified wave refraction computer program was used to construct 
the refraction diagrams used in calculating the simulated spectra. 

This program computes the wave orthogonals (perpendicular to the wave 
crest) based upon the analogy with Snell's law for optics given in 
equation (7)* The program was modified (Poole 1976b) to initiate 
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ray computations at deepwater locations lying along a line perpendicu- 
lar to the initial ray direction. (For this model, deepwater was con- 
sidered when the local water depth was greater than l/U the wave 
length.) All computations were performed at fixed time increments 
referenced to the same initial time. A small family of -wave rays were 
■used, clustered about the initial wave ray and spaced 1 nautical, mile 
apart. The initial ray was chosen to correspond closely with the .data 
flightline which lies along the direction of wave propagation. Based 
on a least squares fit to the inertial navigation data, this angle was 
59° as measured from true North. In the model coordinate system this 
is the angle between the flight track and the X-axis (.Figure 8a). 
Refraction diagrams were computed for wave periods ranging from 6 to 
llj seconds. There was negligible refraction for periods of less than 
6 seconds and little energy associated with periods greater than 14 

t 

seconds, based on examination, of the reference spectrum. Sample 
refraction diagrams ' are presented in figure 8 for wave periods of 6, 

8, 10, 12 and l‘U seconds'. The circular symbols on each figure indicate 
the position of the inertial navigation points. 

These diagrams were then used to compute the energy amplification 
function for each data segment* The normalized energy amplification 
function was computed from equation (ll) for each ray pair that 
hounded the flight path* These values were then averaged over the 
flight path for each data segment* When more than one pair of rays 
hounded the flight path, weighted averages were used* For amplifica- 
tion ^values corresponding to frequencies higher than O.I67 hertz 
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(6 second period) the amplifi cation value vas assumed to he equal to 1 
since no refraction takes place. For frequencies less than 0.071 hertz 
(l4 second period) the amplification function was set at the value 
corresponding to 0.071 hertz. These normalized energy amplification 
functions are presented in figure 9 as a function of frequency for each 
data segment 5 through 10. 

Simulated refracted spectra were computed by multiplying the 
deepwater reference spectrum (figure 6) by the energy amplification 
function according to equation (13):, for each data segment 5 through 
10. The resulting simulated spectra are presented in figure 10. 
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RESULTS AMD DISCUSSION 

In figure 21, the deep water sea state at the time of the experi- 
ment is 1 compared with data from Phillips (1968, page 113, figure 4.8) 
for a fully developed sea and with the deepwater reference spectrum from 
Poole (1976b). The wave field associated with the passage of tropical 
storm Amy contained a significant portion of wave energy at frequencies 
lower them that measured by Poole, and an order of magnitude higher in 
spectral density. The reference spectrum for the Amy study (from figure 
6) is for an area nearly 400 nmi from the center of the storm, but prior 
to the region of shoaling and refraction effects that occur further in- 
shore. The data from Phillips illustrates the equilibrium range of the 
frequency spectrum for wind generated waves. The slope of the spectral 
peak is shown for only three cases in addition to the deepwater reference 
spectra; otherwise only the saturated part of each spectrum is shown. 

For the fully developed sea in figure 11, each possible frequency band 
in the spectrum is present with a maximum amount of spectral energy 
(Pierson, Neumann, and James, 1955 j page 4l ) . The shape of the reference 
spectrum does not match that for the fully developed sea, as would be 
expected since the reference spectrum is not in the generating area. At 
this distance from the storm, the waves have been filtered, with the 
longer waves (lower frequencies) arriving at this location before the 
shorter waves (higher frequencies) that were generated at the same time. 
Thus the reference spectrum with which the inshore spectra are compared 
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is based on a sea state that is not fully developed but is characteristic 
of the filtered sea that is outside the storm area. 

The deepwater reference spectrum, figure 6, is the result of 
averaging the first four deepwater spectra to obtain a greater degree 
of statistical stability than the individual spectrum would have. The 
number of degrees of freedom goes from llU for the individual spectrum 
to 458 for the combined reference spectrum and narrows the width of the 
confidence bands associated with the data. These -bands imply that for 
repeated tests, the spectral density values would lie between 0.9016 
and 1.086 times the values of the reference spectrum 80 percent of the 
time. Changes of the spectral levels that fall within these limits can 
be attributed to random fluctuations in the data, but changes which fall 
outside these values are due to either shoaling or refraction or a combi- 
nation of these effects. Each of these factors effect the energy density 
levels by their effect on wave height according to equation (9)« The 
experimental and simulated spectra are compared with the reference spec- 
trum .in figures 7 and 10, and with each other in figure 12. 

In figure 7> the experimental ' spectra show a gradual decrease in the 

, ^ * * * 

peak spectral density levels to a minimum' in segment 7 (figure 7c) and 
rises again to a maximum in. segment 9 (figure 7e) then drops rapidly to 
its lowest level in the most inshore segment, segment 10 (figure 7f ) • 
Segment 10 is also double peaked. The experimental spectra are further 
characterized by a shift of the frequency at which the peak energy density 
occurs to lower frequency levels as the spectra move shoreward . It is 
in spectral segments 7, 8 and 10 (figures 7c, 7<1, 7f)' that variations in 
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the data occur that exceed the 80 percent confidence hands of the 
reference spectrum and can be attributed to the effects of shoaling 
and/or refraction and not random fluctuations in the data. 

In figure 10 the simulated spectra are compared with the reference 
spectrum. The peak spectral density level shows a general decrease to a 
minimum in spectral segment 8 (figure lOd) , increases to a value greater 
than the reference spectrum in segment 9 (figure lOe) and then has a slight 
decrease in segment 10 (figure 10f). The spectrum becomes double peaked 
in segment 8 (figure lOd) and indicates a shift in the frequency of peak 
spectral density to a lower frequency in segment 10 (figure lOf). Even 
though this fits the general description of the experimental spectra, it 
is only in spectral segment 8 (figure lOd) that the change is sufficiently 
large to fall outside the 80 percent confidence bounds and would predict 
changes which could be considered statistically significant. 

In figure 12, the simulated spectra are compared with the experi- 
mental spectra- Both show the same undulating variation in peak spectral 
values with a minimum near segments 7 and 8, rising again in segment 9 a nd 
decreasing again in segment 10- The overlapping confidence bands in all 
cases except segment 10 indicate good agreement between the simulated and 
experimental spectra for all but the most inshore segment. In figure 13 
the total energies are compared and as in the spectral comparisons , the 
energy levels of the simulated spectra follow the general trend of the 
experimental total energies. A comparison is made of the frequencies at 
which the peak spectral density occurs in figure lU. There is a general 
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decrease in the experimental value from 0-1062 hertz (9*^2 second period) 
in the reference spectrum to 0-0912 hertz (10. 96 second period) in spec- 
tral segment 9- This shifting of the peak wave energies toward the lower 
frequencies (longer wave periods) as the waves move inshore was not 
duplicated by the simulated spectra. With this exception, the theoretical 
technique described in Poole (1976b), when applied to the data from 
tropical storm Amy was able to simulate the changes in the spectra caused 
by shoaling and refraction over most of the shelf area. 

An examination of the behavior of the shoaling and refraction co- 

‘ * * * * * 

efficients over the test region may’ help to understand the observed 

p 

spectral changes in 'the wave field as it moves shoreward. The variation 
of the shoaling coefficient with bottom depth, as shown in Kinsman .(1965 , 
page 159) for the wave frequencies and depths associated with the test 
indicate that the shdaling coefficients for the first 9 segments were 
all less than 1- The effect of this coefficient is to decrease the 
spectral density as the wave field moves inshore, reaching a minimum 
value near spectral segments 8 and 9- Only in segment 10 would the 
shoaling cause an increase in the spectral density levels* Therefore, 
any increase in energy levels in the first 9 segments higher than the 
reference level would have to be due to the effects of the refraction 
coefficient. This value is dependent on the separation of wave rays that 
bracket the data segment and can have values which range from less than 
1 for diverging rays to greater than 1 for converging rays. The energy 
amplification function shown in figure 9 best illustrates the combined 
effect of the shoaling and refraction coefficients in shifting the energy 



density levels as the wave field moves inshore. After an initial 
reduction in values (seen in figures 9 a and 9t>) for the two most off- 
shore segments (segments 5 and 6) there appears to he a general in- 
crease in the amplification function for lower frequency waves as the 
spectral segments move shoreward. This would he as expected since the 
longer waves (i.e., lower frequency waves) are the first to he affected 
hy the bottom and show signs of shoaling and refraction. As the wave 
field moves inshore these effects spread toward the shorter waves (i.e., 
higher frequency waves) which also begin to he affected hy the bottom 
and this shoreward progress can he observed in the amplification function. 
These observations are consistent with the shift of energy to lower 
frequencies seen in figure 7 for the experimental spectra, to a lesser 
degree in figure 10 for the simulated spectra and with the shift of 
the frequency of peak energy density shown in figure Ik. 

The spectra of segment 10 are unique relative to the spectra of the 
other data segments. The experimental spectrum is based on only slightly 
over half the total number of data points that were available for the 
offshore spectra, and about 8 percent of these were taken through the 
surf zone. The assumption of homogeneous wave conditions is also pro- 
bably violated to some degree in this region as this is the area of 
maximum changes in the waves due to shoaling and refraction effects. 

For these reasons , the resulting spectra might be expected to show some 
anomolies. Yet the double peaked nature of the experimental spectrum 
is similar to the spectra recorded by the U.S. Army Corp of Engineers 
Research Center f s (CERC) .shoreline gauge located at Jennette's Pier, 
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Nags Head, North Carolina (Harris, 1975 ), and shown in figure 15- The 
experimental flight was coordinated with CERC personnel and was designed 
to pass directly over the recording station. The CERC spectra were 
generated from data recorded 'before, during, and after the flight passed 
overhead. Spectra 15b and 15c bracket the overflight. Nearly all the 
spectra have similar double peaks with the major peak near 0.08 hertz 
and the secondary peak near 0.15 hertz. This compares with 0.085 hertz 
and 0.153 hertz, respectively, for the experimental spectrum. As an 
indication of total energy, a significant wave height of 1 . 58 m was com- 
puted from the experimental spectrum by trapezoidal integration of the 
area under the curve (under the assumption of a Guassian distribution 
of wave height). This compares well with the 1.86m and l.jbm for the 
CERC spectra nearest to flight time. Even though measured by two 
different techniques, the experimental spectra and the CERC spectra 
have similar- characteristics . This would appear then to be an accurate 
description of the energy distribution in the most inshore area. 

The secondary peaks observed in the segment 10 spectra, and not the 
offshore spectra, appear to be caused by nonlinear effects. In the 
rapidly shoaling waters just off the beach, the distortion of the wave 
form and the increase in wave height due to convergence creates nonlinear 
changes in the phase speed which are not of as large a magnitude in the 
deeper offshore waters (Grosch and Comery, 1977)- In nonlinear wave 
theory these changes increase the terms which are direct multiples of the 
wave frequency (Grosch, 1978). In the experimental spectrum of segment 
10 , the secondary peak is at a frequency nearly double that of the primary 
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peak frequency. In the CERC spectra shown in figure 151 , peaks can he 
seen which correspond to multiples of 2, 3, and h times the primary 
peak frequency. 

Although the experimental spectra from the aircraft data compared 
well with the CERC spectra, the simulated spectra failed to match this 
distribution (figure 12f). The significant wave height of 1.79m compares 
well with 1 . 70 m for the experimental spectrum, but the simulated spectrum 
is single peaked, and at a much higher level of energy density. In this 
region of high wave activity, a number of assumptions on which the simu- 
lated spectra are based begin to fail. The small amplitude assumption 
requires wave height to be small relative to wave length and water depth 
(Kinsman 1965 , pl25) but in the shoaling water the wave lengths shorten 
and wave height increases: Also the equations' used to compute refraction 

are strictly valid only for a constant water depth but have been used 
sue cess fully "in - areas with a mild slope (Dobson, 1967 )* In segment 10, 
since the depth goes to 0, the slope can no longer be considered mild. 

For these reasons , the simulated spectrum might not be expected to match 
the experimental results. 

It is an assumption of wave refraction theory that the energy between 
wave rays remains constant and, does not cross the ray boundaries . Thus 
when rays diverge, the energy density decreases, and when they converge, 
the density increases. If the wave rays cross, a "caustic" is formed, 
indicating that the waves in this region have become so high as to be 
unstable and are breaking (Kinsman 1965 , p 158). The practical applica- 
tion of this theory requires a decision on how to treat the energy 
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associated with the wave rays which formed the "caustic”. In most of 
the offshore spectral segments , few rays crossed, and the energy from 
these rays was included in the calculation. The spectra in data seg- 
ment 10 (figure lOf and 12f) were also computed in this manner. However, 
since in segment 10 there were many "caustics”, this spectrum was also 
computed "by a second method in which the energy from the crossed rays 
was not included in the computations. The results of this method are 
shown in figure 1 6 . The simulated spectrum using this technique "becomes 
double peaked and provided a better match with the experimental spectra 
than did the standard technique. For the same reasons that figure 12f 
did not match the experimental results, figure 1 6 should not either, but 
did. Whether this is a technique that better simulates the loss of 
energy due to wave breaking in this region, or simply a fortuitous 
result cannot be assessed from this single sample. That assessment will 

await repeated tests and is presented here as a point of interest for 

♦ 

* 

future work in this area. 

Had the wave pattern displayed a more clearly defined direction of 
swell propogation or had there been only a single well defined swell, 
the overall comparison between the simulated and experimental spectra 
would possibly have been improved. Such a condition would better match 
the monochromatic conditions assumed in the theoretical model, and the 
assumption of wave movement along the flight path, than did the conditions 
measured. Photographs from the flight (Morris 1978) show a number of wave 
patterns propogating from the general direction of the storm. Even 
though there appears to be a slight dominance of the off flightline waves. 



28 


this could he attributed to the sun angle highlighting the photographs. 

The initial flight angle was chosen by visual observation from the air- 
craft at flight time, and even though a wave pattern similar to that seen 
in the post flight photographs was observed, the dominant direction used 
in the refraction model was the direction chosen for the flight path. 

The local surface winds which were present were approximately per- 
pendicular to the flightline and would not have greatly influenced the 
wave measurements along the direction of flight. 

The amplification function curve for each spectral segment was 
constructed from over 20 sets of refraction diagrams to ensure that each 
major fluctuation in the curve was included. This is more detail than the 
modeling technique requires. The model computes the values of simulated 
spectral density for the frequencies which correspond to those from the 
reference spectrum (such as those shown by the arrows located along the 
abscissa of figure-9d). It is clear that much activity can occur between 
these points and will not be treated by the model, especially at the 
lower frequencies. As long as the computational points follow the general 
trend of the amplification function, the results would appear to be re- 
presentative. In cases where the energy spectrum is spread over a broad 
frequency range, and where lower total energy is reflected in a less 
active amplification function, the modeling technique would not be overly 
sensitive to the location of the computational points. As energy concen- 
trates toward the lower frequency range, however, the model will attempt 
to define the increasing concentrations of enehgy with a decreasing 
number of points. To increase the number of frequencies used in this 
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model* the digitization rate will have to be increased as this will 
better define both the reference spectrum and the simulated spectra. 

It should be noted that this can cause large increases in the model* s 
computer requirements. 

Although the data sampled does not represent the ideal conditions 
required for verification of the theoretical modeling technique* the 
wave climate is probably typical of the conditions for which it can be 
used. The study by Poole (19?6b)* illustrated in figure II* apparently 
represents a minimum sea condition for which the model can be applied. 

The present study, with more energy at the lower frequency has estab- 
lished a level at which the model shows a sensitivity to the shoaling 
and refraction effects which are recognizable from the random fluctuations 
of wave energies. Studies of wave conditions which fall within the 
bounds established by the fully developed sea and the present study 
should help further- establish the degree to which the model is sensi- 


tive to these changes. 
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CONCLUDING REMARKS 

This report presents a comparison of simulated wave spectra with 
experimental wave spectra. Experimental spectra were computed from 
surface elevation measurements taken from an aircraft over the continental 
shelf region off North Carolina after passage of tropical storm Amy on 
July 1975. Simulated spectra corresponding to 6 of the experimental 
spectral data segments were computed using a modified linear wave 
refraction model. 


Conclusions 

The deep water wave field associated with tropical storm Amy was 
not that of a fully developed sea, hut probably typical of a wave field 
outside the wave generating area. High frequency waves had been filtered 
and energy concentrated toward the low frequencies. When compared with 
the spectrum of “the deepwater wave field, the inshore experimental 
spectra showed a statistically significant change in spectral levels for 
several of the inshore spectra which can he attributed to shoaling and 
refraction effects. The spectra also showed a consistent shift of energy 
toward the lower frequencies as the wave field moved shoreward. The most 
inshore spectra agreed with the results of CERC shoreline spectra 
recorded simultaneous to the experiment. The simulated spectra followed 
the shift in energy levels of the experimental spectra for all hut the 
most inshore spectral segment. The total energy levels were within 20 
percent of the experimental levels hut the shift in frequency of peak 
spectral density was not predicted hy the simulated spectra. 
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Given a deepwater spectrum, the theoretical techniques used to 
simulate the spectra over the shelf appears to he capable of predicting 
the changes in spectral shape and energy levels that occur due to shoal- 
ing and refraction effects to within 5 nmi of the shoreline. Shoreward 
of this region the simulation technique fails to match the energy distri- 
bution of the experimental spectrum. Future applications of this tech- 
nique for higher energy levels, will require careful selection of the 
rate at which the data is digitized to ensure a sufficient definition 
of the experimental spectra and to ensure the ability of the theoretical 
modeling technique to simulate the changes due to shoaling and refraction 
effects. 


Suggestions for Future Work 

1. Storm systems stronger than Amy will produce wave fields of even 
more interest to those concerned with wave energy distributions in 
the inshore area. The results presented in this report indicate the 
ability of the model to simulate the inshore spectra for a moderately 
strong storm such as Amy, but results will be needed from systems con- 
centrating even higher energy levels to help determine the sensitivity 
of the modeling technique. The spectral comparison with fully developed 
sea spectra (figure 11) establishes bounds that can be used in evaluating 
the value of using a deepwater spectrum for simulation. 

2. An assumption of the technique used to determine the experimental 
spectra is that the data is gathered along the direction of swell pro- 
pogation. This is difficult to determine while observing the sea from an 
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aircraft and forecast wind information is often used to determine the 
flight direction. A comparison of spectra from data collected from 
several flight directions for the same wave field would help determine 
the sensitivity of the technique to errors in choosing the direction 
of swell propagation. 

3. One impediment to the use of this theoretical technique to simulate 
the spectra is in determining the amplification factors for each spectral 
segment. The present system requires the use of an overlay to determine 
which ray pairs bracket the flight path and the manual summing and averag- 
ing of the amplification functions associated with these rays. This is 

a most tedious job and an obstacle to rapid application of this technique. 
The procedure should lend itself to at least partial automation which 
would greatly enhance the application of this technique and promote its 
use. 

4. Although the - method for simulating spectra predicts the energy dis- 
tributions across the shelf, it fails in the most inshore case where 
there is frequencly even more need for the results. On the assumption 
that the wave conditions just offshore are indicative of the conditions 
at the beach, it should be of interest to see just how near to the surf 
zone the spectral segment can be extended and still simulate the experi- 
mental spectra. In the present case the simulation matches to within 
4.7nmi of the beach. By selectively breaking the spectral segment to 
form a moving spectral band it should be possible to determine how far 
inshore the segment could be extended before failing. The nearer to 
shore the spectra matched, the higher would be the confidence that the 
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results would indicate conditions along the beach. 

5* The technique described in this report, which excludes from the 
most inshore spectra the energy of those waves which crossed and formed 
a "caustic", should be explored. The uniqueness of this result can only 
be verified by application to another independent set of data. If the 
result is other than chance, the techniques would offer a method of 
extending the simulated spectra close to the shoreline. 
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APPENDIX 

SYMBOLS 

b ray separation distance, meters 

C autocorrelation function 

c celerity, meters /second 

d water depth, meters 

2 

E energy density, meters 

F normalized energy amplification function 

f circular frequency, hertz 

2 

g gravitational constant, meters/second 

H wave height, meters 

K refraction coefficient 

-r 

K g shoaling coefficient 

k wave number 

& counter (0 < Ji < M) 

L wave length, meters 

M total number of lags> 

N total number of sample data points 

2 

S spectral density, meters /hertz 

t time, seconds 

V velocity, meters/second 

a angle formed by wave crest and local depth contour, degrees 
r) wave elevations profile, meter 

radian frequency, radians/second 


0) 
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T number of lags 

Subscripts : 
a apparent 

0 initial (deepwater) value 

A bar over a symbol indicates an average* 
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TABLE 1 

GROUND TRACK OF TROPICAL STORM AMY 


DATE 

TIME 

LATITUDE 

LONGITUDE 

MAX. WIND STAGE 

June 27 

OOOOZ 

27. 5N 

79- ow 

25 KT Trop. Dep. 


0600Z 

28.5 

79-0 



1200Z 

29.5 

79-0 



1800Z 

30.5 

79-0 


June 28 

OOOOZ 

31.5 

78.8 



0600Z 

32.1; 

78.7 



1200Z 

33.3 

78.0 



1800Z 

34.0 

77.0 

30 

June 29 

OOOOZ 

34.4 

75.8 

35 Trop . Storm 


0600Z 

3b. 0 

74.8 

40 


1200Z 

33-8 

73.8 

45 


1800Z 

33.8 

72.8 

50 

June 30 

OOOOZ 

34.3 

71.6 



0600Z 

35.6 

70.8 

55 


1200Z 

35-9 

70.5 

60 


1800Z 

36.2 

70.2 


July 1 

OOOOZ 

36.2 

69*8 



0600Z 

36.2 

69.4 



1200Z 

36.2 

68.3 



1800Z 

38.7 

67.2 


July 2 

OOOOZ 

37. 4 

66.7 



0600Z 

37.3 

65.9 



1200Z 

37.3 

65.1 



1800Z 

37.3 

64.1 


July 3, 

OOOOZ 

37-7 

62.8 

55 


06OOZ 

38.2 

61.2 



1200Z 

39-3 

59-6 



1800Z 

40.5 

58.0 

50 
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TABLE 1 (Continued) 


DATE 

TIME 

LATITUDE 

LONGITUDE 

MAX. WIND 

July 4 

ooooz 

42-5 

54. 8 



0600Z 

44.5 

51.6 



1200Z 

1*7-0 

48.0 

45 


STAGE. 


Extra- 

Tropical 
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TABLE ,2 

FLIGHT NAVIGATION LOG 

TIME (EDT) LATITUDE LONGITUDE 


Hr. 

Min. 

Sec. 

Degrees' 

Minutes 

Degrees 

Minutes 

16 

12 

00 

36 

33. 9N 

74 

15 . 4w 

1 6 

l4 

00 

- 

- 

- 

- 

16 

16 

00 

36 

28.8 

7b 

27.2 

16 

18 

00 

36 

26.3 

7b 

32.7 

l 6 

20 

00 

36 

23.3 

7b 

38.7 

16 

22 

00 

36 

20 . 4 

7b 

1 * 1*. 6 

16 

24 

00 

36 

IT. 3 

7b 

51.1 

16 

26 

00 

36, 

,l4.8 

7b 

56.5 

16 

28 

00 

36 

11.3 

75 

04.1 

16 

30 

00 

36 

09.3 

75 

08.0 

16 

32 

00 

36 

06.4 

75 

13.7 

16 

34 

00 

36 

03.5 

75 

14.7 

16 

36 

00 

36 

00.9 

75 

25-0 

16 

39 

26 

35 

55.1 

'75 

34.1 
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Figure 1.- Ground track of tropical storm Amy. 




Figure 2. - Photograph of tropical storm Amy, 1500 GMT, July 2, 1975, 
taken from the Synchronous Meterological Satellite-1. 
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• (a) Segment 1 

Figure 5. -.Experimental deep-water spectra for flight track segments. 
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Figure 5.- Continued 
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Figure 5.- Concluded 
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Frequency, Hz 

Figure 6.- Deep-water reference spectrum. 
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Frequency, Hz 
(a) Segment 5 

Figure 7.- Comparison of the experimental spectra with the reference spectra 
for each data segment. 
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(e) Segment 9 
Figure 7.- Continued 
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(b) Segment 6 
Figure 10.- Continued 
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(c) Segment 7 
Figure 10.- Continued 




Spectral density, M**2/Hz 


71 



Frequency, Hz 
(d) Segment 8 


Figure 10.- Continued 
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Frequency, Hz 
(e) Segment 9 
Figure 10.- Continued 
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Figure 11.- Comparison of reference spectra with the frequency 
spectrum of wind generated waves in the equilibrium 
range. From Phillips (1968, page 113, figure 4.8). 
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Frequency, Hz 
(a) Segment 5 

Figure 12.- Comparison of experimental and simulated wave spectra for each 
.data segment. 
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Frequency, Hz 
(e) Segment 9 
Figure 12.- Continued 
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(f) Segment 10 
Figure 12.- Concluded 
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Spectrum number 

Figure 14.- Experimental- and simulated frequency of peak spectral 
density variation with data segment. 
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Figure 15.- Wave spectra from the CERC shore gauge at Jeanette* s Pier, Nags Head, 
North Carolina (Harris, 1975). Significant wave heights are shown 
for each spectra. 
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